Universal probes for antiferromagnetic correlations and entropy 
in cold fermions on optical lattices 



E. V. Gorelik, 1 D. Rost, 1 T. Paiva, 2 R. Scalettar, 3 A. KKimper, 4 and N. Bliimer 1 

1 Institute of Physics, Johannes Gutenberg University, Mainz, Germany 
2 Instituto de Fisica, Universidade Federal do Rio de Janeiro, Brazil 
3 Department of Physics, UC Davis, USA 
4 University of Wuppertal, Wuppertal, Germany 
(Dated: June 6, 2012) 

We determine antiferromagnetic (AF) signatures in the half-filled Hubbard model at strong cou- 
pling on a cubic lattice and in lower dimensions. Upon cooling, the transition from the charge- 
excitation regime to the AF Heisenberg regime is signaled by a universal minimum of the double 
occupancy at entropy s = S/(NkB) — s* ~ ln(2) per particle and a linear increase of the next-nearest 
neighbor (NNN) spin correlation function for s < s* . This crossover, driven by a gain in kinetic 
exchange energy, appears as the essential AF physics relevant for current cold-atom experiments. 
The onset of long-range AF order (at low s on cubic lattices) is hardly visible in nearest-neighbor 
spin correlations versus s, but could be detected in spin correlations at or beyond NNN distances. 
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Materials with strong electronic correlations are, due 
to their increasing technological importance, a prime 
subject of current research [l|, 0|. Theoretical inves- 
tigations of corresponding Hubbard-type models have 
shed light on many strong-coupling phenomena, includ- 
ing metal-insulator transitions, non-Fermi-liquid behav- 
ior, and various types of magnetic and orbital order Q. 
However, important questions remain open, most notably 
regarding high-temperature superconductivity, for which 
so far no mechanism could be established. Recently, a 
novel class of correlated Fermi systems, namely ultracold 
fermionic atoms (such as 40 K and 6 Li) on optical lat- 
tices, has opened a new promising direction of research: 
Cold atoms are predicted to serve as quantum simulators 
for the Hubbard type solid-state Hamiltonians of interest 



Indeed, the Mott metal-insulator transition has re- 
cently been demonstrated in two-flavor mixtures of 40 K 
on cubic optical lattices by experimental observation and 
quantitative theoretical analysis of signatures in the com- 
pressibility [3] and in the double occupancy Q. This 
success established that the single-band Hubbard model 
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(with hopping amplitude t, on-site interaction U , and 
hi a = c\ a c i(J ) can be realized to a reasonable accuracy 
using ultracold fermions in the interesting interaction 
range, which certainly supports the hopes of accessing 
also less understood Hubbard physics in similar ways. 

However, attempts to realize and detect quantum mag- 
netism in cold lattice fermions have proven extremely 
challenging. In fact, it is even difficult to verify spe- 
cific signatures of antiferromagnetic (AF) correlations 
which are believed to play an important role in high- 
temperature superconductivity. This type of physics 
clearly has to be under control before cold fermions on 



optical lattices can really play a useful role as quantum 
simulators of materials with strong electronic correla- 
tions. The failure to detect AF signals has primarily 
been attributed to cooling issues Indeed, the cold- 

est experiments for repulsive fermions on optical lattices 
have thus far reached central entropies per particle of 
s = S/(Nk B ) « ln(2) w 0.7 [HG3 while AF long-range 
order (LRO) on an isotropic cubic lattice is expected only 
for entropies s < s N w 0.4 0, E3] - 

An important feature of cold-atom systems is their in- 
homogeneity, induced by the trapping potential. On the 
one hand, this is beneficial for quantum magnetism, since 
entropy is effectively pushed out of a half-filled core: this 
aspect is a major theme of current research [ll|, ■ 
On the other hand, any AF region will necessarily be of 
limited spatial extent (l8| . so that the thermodynamic 
concept of LRO is not fully applicable. In addition, as 
we will show, the nearest-neighbor (NN) spin correlation 
function, which is currently addressed (using modulation 
spectroscopy (qL fl9t or super lattices 0) in AF related 
experiments, is hardly sensitive to the onset of LRO even 
in the thermodynamic limit. 

In this situation, one may ask the following: (i) Is 
there a threshold distance beyond which spin correlations 
have "long-range characteristics" and (ii) can we define 
"finite-range antifcrromagnctism" as a unique scenario 
with universal properties, appearing only in a certain en- 
tropy range? The answer to both questions is "yes" : The 
essential AF correlation physics emerges already at en- 
tropies s < m(2), i.e. in reach of current cooling tech- 
niques. Since, in addition, the threshold distance turns 
out to be rather small (but larger than one lattice spac- 
ing), our quantitative predictions should enable experi- 
mentalists to verify specific AF signatures with current 
system sizes, i.e., to get the long-sought grip on quantum 
magnetism. 

In the following, we discuss first an enhancement of the 
double occupancy D (i.e. also of the interaction energy 
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Ei ut = D U) at low temperatures T which has previously 
been proposed as an AF signature on the basis of dy- 
namical mean-field theory (DMFT) DMFT results 
for a half-filled cubic lattice at strong coupling U/t = 15 
are compared with direct determinantal quantum Monte 
Carlo (DQMC) [12, simulations. The comparison is 
extended to the full dimensional range based on DQMC 
and Bethe ansatz (BA) [24| data in dimensions d = 2 
and d = 1, respectively. High-precision estimates of the 
entropy s(T) allow us to switch to the experimentally 
relevant entropy representation. An asymptotic collapse 
of the curves D(s) is observed with respect to dimen- 
sionality, with universal minima at s* f=a ln(2), and no 
significant features at sn in the cubic case. Additional 
specific signatures of finite-range AF order are found in 
the kinetic energy and in spin correlation functions, with 
different degrees of universality. Finally, the perspectives 
for detecting LRO are discussed using stochastic series 
expansion (SSE) results for the Heisenberg model. 

AF signatures in the double occupancy. Accord- 
ing to DMFT, the low-T formation of an AF core in a 
fermionic cloud on an optical lattice (with central half fill- 
ing, n = 1) is signaled, at strong coupling, by a distinct 
enhancement of D in the same region [2l| . As a func- 
tion of temperature, DMFT predicts nearly fiat curves 
D(T) in the range T > T$ MFT , i.e. above its estimate of 
the Neel temperature, and a sharp increase below, with 
a kink and absolute minimum at T® MFT . This is clearly 
seen, for U/t = 15, in Fig. [lja) (circles). The absolute 
low-T increase of D is largest for U/t rj 12; it should be 
detectable, according to real-space DMFT, even in exper- 
iments integrating over the inhomogeneous cloud [21] . 

Not all aspects of this DMFT scenario are, however, 
realistic: After all, DMFT is exact only in the limit of 
infinite coordination number Z — > oo (with Z = 2d for 
hypcrcubic lattices) and overestimates the Neel tempera- 
ture by up to 30% in the simple cubic case dHU. Thus, 
the sharp kink in D{T) seen in Fig-CfTa) at T r / MFT sa OAt 
cannot be physical. One might expect a shift of the 
DMFT results toward lower temperatures, as well as 
some broadening in the cubic case and more radical 
changes (at least) for d < 2; only at high temperatures 
does the accuracy of DMFT estimates for D follow al- 
ready from series expansions (in d = 3) [27I . 

Impact of dimensionality. - Indeed, DQMC estimates 
of D(T) [diamonds in Fig. QTa)] agree with DMFT for 
T/t > 1 within error bars, which are smaller than sym- 
bol sizes [2q . Surprisingly good agreement is observed 
also at T/t < 0.3. As expected [Hi, the DMFT kink 
is smeared out in the DQMC data toward a broad mini- 
mum. At suitably rescaled [HJ interactions, DQMC data 
for d = 2 (squares) show remarkably similar behavior. 
Only the case d = 1 (dotted line) deviates more drasti- 
cally at intermediate and low T. Note that the position 
of the minimum in D(T) shifts upward with decreasing 
d, i.e., opposite to the naive expectation. 

As optical lattice and interactions are switched on for 
the ultracold atoms in a nearly adiabatic process [ll[ , the 
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FIG. 1. (Color online) Hypercubic lattice (1 < d < 3) at 
strong coupling U/(\f~Zt) ~ 6: (a) Double occupancy D(T) 
as estimated from DMFT (d = 3, circles), DQMC (d = 3: 
diamonds; d = 2: squares), and BA (d — 1, dotted line), (b) 
Corresponding estimates of entropy per particle s — S/(NkB). 



entropy s (and not T) is the experimentally relevant con- 
trol parameter. Fig. [TJj shows numerically exact data for 
s(T), obtained directly for d = 1 and via the thermody- 
namic relation £(/3) = ln(4) + /3E(p) - f£ dj3'E(P') [with 
;3 = l/(k B T) and energy E] for d = 2,3 and DMFT. 
Again, the agreement between d = 2 and d = 3 is strik- 
ing; the latter results converge to the Heisenberg limit 
for T < 0.8t. Remarkably, the proper DMFT solution 
(circles) is close to the DQMC result for the cubic lat- 
tice (diamonds) at T < 0.3t; only the metastable non- 
magnetic DMFT solution (thin solid lines), considered in 
previous studies [27| . remains far off. 

Figure [21 obtained by combining the data of both pan- 
els of Fig. [TJ conveys our first central message: As a func- 
tion of entropy, the double occupancy is surprisingly uni- 
versal at strong coupling, with a minimum at s* ~ ln(2) 
in all dimensions and generally similar shapes. At con- 
stant rescaled interaction U, the curvature around the 
minimum increases with increasing dimensionality un- 
til it becomes sharp in the DMFT limit, where it cor- 
responds to the Neel transition. As seen in the inset (for 
d = 1), the minimum becomes also sharp and approaches 
s* at constant dimensionality in the strong-coupling limit 
U — > oo. Evidently, the minimum in D(s) separates two 
regimes with quite different physical properties: (i) the 
regime s > s*, which smoothly approaches the Hartree 
limit (with D = (n t )(rij_) = 0.25) for T -> oo, and 
(ii) a low-temperature regime, with no discernible sub- 
structure. It is clear that the latter regime must be char- 
acterized by spin coherence, since s < ln(2) can occur in 
a two-flavor system at n = 1 only by the development of 
(possibly short ranged) magnetic correlations. 

In fact, any positive deviation of D(s) from the non- 
magnetic background (shaded in Fig. [2]) should be linked 
to AF correlations, generalizing Takahashi's ground state 
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FIG. 2. (Color online) Hypercubic lattice at strong coupling: 
Double occupancy vs entropy. In all cases, the minimum of 
the double occupancy corresponds to s* w ln(2) (dotted line). 
The shaded area indicates the nonmagnetic contribution to D. 
Inset: D(s) in d = 1 for various interactions. 
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Here, {&i ■ <Tj)o 618 is the nearest-neighbor correlation in 
the quantum Heisenberg model at T = (for Pauli ma- 
trices <t), which is stronger in lower d: 
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-1.77 (d = 1) [34], 
-1.34 (d = 2) [35], 



-1.20 (d = 3) [36], 
-1.00 (d = oo), 



consistent with our finite- T results. Thus, irrespective 
of the measurement technique, signatures of AF corre- 
lations may be easier to detect experimentally, at fixed 
s, for lower (effective) dimensionality. Conversely, a tun- 
ing of the hopping amplitude in z direction could help to 
discriminate magnetic effects from those of charge excita- 
tions; similar ideas involving frustration will be explored 
in a separate publication [371 ] . 

Energetics and spin correlations. - Up to a numerical 
pref actor (of 15/ v6), the curves in Fig. [2] represent the 
rescaled interaction energy E- m t / '(y/Zt) . Complementary 
AF signatures appear in the kinetic energy, shown in Fig. 
J3ja). In particular, Ekin(s) exhibits a maximum at s s* 
for d > 3 [while thermodynamic consistency ensures a 
monotonous total energy E(s) = Ei nt (s) + E^i D (s)}. The 
associated redistribution of quasi momentum at s < s* 
might be a worthwhile experimental target. 

As mentioned in the Introduction, nearest-neighbor 
(NN) spin correlations are, so far, playing a central role in 
the experimental quest for AF signatures in cold lattice 
fcrmions. However, as seen in the lower part of Fig.[3jb), 
the NN spin correlation functions (symbols and dashed 
and dotted lines) have strong high-entropy tails in all 
physical dimensions d < 3, with only trivial features at 
s w s* , below which the Heisenberg model (solid lines) 
becomes applicable. Conclusions concerning the presence 
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FIG. 3. (Color online) Hypercubic lattice at strong coupling: 
(a) Rescaled total and kinetic energy vs entropy, (b) Spin 
correlations (<Xj ■ &i+s) between nearest [8d=i = 1, <5<i=2 = 

(1.0) , <5d=3 = (1,0,0))] and next-nearest [6d=i = 2, Sd=2 = 

(1.1) , <5 d=3 = (1,1,0)] neighbors. Inset: NNN and infinite- 
range correlations for Heisenberg model (d = 1 from Ref. 



or proximity of AF order could be drawn from corre- 
sponding experimental data only via theoretical look-up 
tables. 

In contrast, the next- nearest- neighbor (NNN) spin cor- 
relation functions are essentially zero for s > s* and take 
off linearly below in all dimensions. So the mere presence 
of significant NNN correlations already implies that the 
Heisenberg regime s < ln(2) has been reached. Remark- 
ably, the results for d — 1 and d — 2 arc indistinguish- 
able (and close to DMFT) for s > 0.5; only those for 
d = 3 arc slightly below. The same picture emerges in 
the Heisenberg model (upper set of curves in the inset of 
Fig. [3|), with slightly larger absolute values. This near- 
universality with respect to d establishes that the cross- 
ing point between short-range physics, where spin corre- 
lations increase with lowering d, and long-range physics, 
where spin correlations decay quickly toward low d (and 
remain finite in the limit of distance 8 — > oo at s > 
only for d > 3 [3!| , cf . the lower set of curves in the inset 
of Fig. [3]) is essentially at the NNN distance. 

In this sense, spin correlations at or beyond the NNN 
distance are much more representative of genuine AF 
physics than the unspecific NN correlations. This is 
true even regarding LRO, as illustrated in the Heisen- 
berg limit (chosen due to the higher achievable accuracy) 
in Fig. 01 While the entropy derivatives of the NN spin 
correlations in d — 3 (squares) are featureless [as those 
of the spin correlations in d = 2 (thin lines)], a distinct 
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FIG. 4. (Color online) Derivatives of NN (squares) and NNN 
(circles) spin correlation functions of the Heisenberg model 
with respect to entropy (main panel) and temperature (inset) 
in d — 3. Thin solid and dashed lines: Heisenberg results in 
d = 2. 

peak at the critical entropy sn for AF LRO in d = 3 [l3[ 
is visible in the NNN data (circles), with a nearly flat 
plateau above. As a comparison with corresponding T 
derivatives (with clear peaks at both in the NN and 
NNN data - see the inset of Fig. Q| shows, this qualita- 
tive difference is associated with the use of the entropy 
as a control parameter. With ex per iments being confined 
to this entropy parametrization [4l| , it is clear that a de- 
tection of the Neel transition via spin correlations would 
require measurements at least at NNN distances. 



Conclusion. - We have disentangled generic aspects of 
antiferromagnetism both from those specific to infinite- 
range order (which is not attainable, by definition, in fi- 
nite atomic clouds) and from trivial nearest-neighbor cor- 
relations that persist even at high temperatures, where 
spin models are not adequate. Our results establish 
that the regime s < s* « ln(2) is characterized by 
"finite-range antiferromagnetism" with remarkably uni- 
versal properties. It may be detected experimentally, at 
strong coupling, by a negative slope in D(s) [or D(T)\ 
or by the onset of longer-range (beyond NN) spin corre- 
lations [HJ. Long-range order appears inessential in the 
cold-atom context and largely decoupled from the basic 
correlation mechanisms (cf. recent iron layer experiments 
[43|): The mean-field critical entropy s* is experimentally 
more relevant than the critical entropy of an infinite sys- 
tem (in d = 3) HU]. Thus, reaching s < sn ~ 0.4 in cu- 
bic systems would not guarantee additional insight and 
global equilibration (the time scales for which may exceed 
experimental life times j4f|) is not essential. Tuning the 
dimensionality (and/or adding frustration) for discrimi- 
nating AF effects appears much more promising. 
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